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ABSTRACT 

The VLA-COSMOS Large Project has imaged the 2 deg^ COSMOS field 
with a resolution of 1.5 arcsec and a sensitivity of about 11 /xJy (la) yielding to 
a catalog of ~ 3600 radio sources. In this paper we present a further analysis of 
the VLA-COSMOS Large Project catalog of radio sources aimed to: 1) quantify 
and correct for the effect of bandwidth smearing in the catalog, 2) determine the 
incompleteness produced by the noise bias and the resolution bias in the new 
catalog and 3) derive the radio source counts at 1.4 GHz. 

The effect of bandwidth smearing on the radio sources in the catalog was 
quantified comparing the peak and total fiux densities in the final mosaic and in 
each of the individual pointings where the source was closest to the center of the 
field. We find that the peak fiux densities in the original VLA-COSMOS Large 
Project catalog have to be divided by a factor about 0.8 or 0.9, depending on the 
distance from the mosaic center. 

The completeness of the radio catalog has been tested using samples of sim- 
ulated radio sources with different angular size distributions. These simulated 
sources have been added to the radio image and recovered using the same tech- 
niques used to produce the radio catalog. The fraction of missed sources as a 
function of the total flux density is a direct measure of the incompleteness. 

Finally, we derived the radio source counts down to 60 /zJy with unprece- 
dented good statistics. Comparison to the flndings of other surveys shows good 
agreement in the flux density range 0.06-1 mJy conflrming the upturn at ~ 0.5 
mJy and a possible decline of the source counts below ~ 0.1 mJy. 
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Introduction 



In recent years, we have experienced a renaissance of deep radio surveys, usually associ- 
ated with multi- waveband observational programmes mainly designed to probe the formation 
and evolution of galaxies and Super Massive Black Holes (SMBH). These surveys, mostly 
carried out at 1.4 GHz, have enabled the study of sub-mJy and /xJy radio source popula- 
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While it is clear that the /iJy population is a mixture of different classes of objects (radio 
loud and radio quiet AGN, starburst galaxies, spirals), the exact contribution of these differ- 



ent classes to the radio source counts is still not ver y well estab 



on the optical lirait of the spectroscopic follow-up (IBenn et al 
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Simpson et al.ll2006l ). The role played by the cosmological e voluti on of the different 



classes of objects is even more uncertain. Recently, ISmolcic et al.l (120061 ) have developed a 
new classification method to separate star forming galaxies and AGN based on photometric 
rest-frame colors for local g alaxy samples. The a pplication of this classification method to 
the VLA-COSMOS sources (ISmolcic et al.ll2008al ) shows that the radio population above ~ 
50 yuJy is a mixture of roughly 30-40% of star forming galaxies and 50-60% of AGN galaxies, 
with a minor contribution (~ 10%) of QSOs. 



The international COSMOS (Cosmic Evolution) survey (IScoville et al.l l2007al ) is de- 



signed to probe the correlated evolution of galaxies, star formation, AGN up to high redshift. 



the X-rays to the radio domain ( 
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2004. 


2007 


) supported by extensive optical spectroscopic campaigns using mainly the VLT/VIMOS 



and the Magellan/IMACS instruments fiLillv et al.ll2007l : iTrump et al.l 120071 ) 



The radio observations of the COSMOS field were done with the VLA in A and C 
configuration and are described in Schinnerer et al. (2004, 2007). The whole 2 deg^ field 
was imaged in mosaic mode with a resolution of 1.5 arcsec. The sensitivity is uniform 
in the inner 1 deg^ with an average rms value of 10.5 /iJy/beam and it increases almost 
regularly to the outer regions of the 2 deg^ field. From the mosaic a catalog of radio sources, 
the VLA-COSMOS Large Project Catalog, was extracted. The catalog lists ~ 3600 radio 
sources selected above the local 4.5cr threshold over an area of 2 deg^. A full description of 
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the observations, data reduction and catalog extraction is given in Schinnerer et al. (2007, 
hereafter SOT). 

In this paper we present a further analysis on the VLA-COSMOS Large Project Catalog 
aimed to quantify the impact of bandwidth smearing (Section 2) and the survey completeness 
(Section 3). Finally, we derive the radio source counts and compare them with those derived 
from similar surveys (Section 4). 

A more detailed interpretation of the radio data is left to other papers which are cur- 
rently in preparation. I n particular, the op tical identifications of the VLA-COSMOS radio 
sources are presented by lCiliegi et al.l (120081 ) . the fraction of star-for ming galaxies and AGN 
in the sub-mJy population of the VLA-COSMOS field is discussed in lSmolcic et al.l (j2008al ) 
an d the dust un-biased c osmic star formation history is derived from the 1.4 GHz radio data 



m 



Smolcic et al.l fl2008b[ ) 



2. Bandwidth smearing correction 

Bandwidth smearing (or chromatic aberration) affects all synthesis observations made 
with a finite width of the receiver channels. Imaging sources at large distances from the phase 
center can yield to radial smearing, increasing the source size and reducing the peak flux 
density, while the integrated flux density is conserved. The image smearing is proportional to 
the fractional bandwidth (the bandwidth divided by the central frequency o f observation) an d 



to the distance from the phase center in units of the synthesized beam (e.g. lThompsonlll999l ) . 
For single pointing observation the bandwidth smearing simply increases with the radial 
distance from the phase centre. For a mosaic image, produced by multi-pointing observations, 
the contribution from all the pointings has to be taken into account and the resulting radial 
pattern will be much more complicated depending also on the adopted spacing pattern of 
the individual pointings. 

In order to quantify the effect of bandwidth smearing on the VLA-COSMOS observa- 
tions we ran the same procedure that produced the radio catalog on each of the 23 individual 
pointings (see S07). For the strongest sources we compared their peak and total flux densi- 
ties in the final mosaic with the corresponding peak and total flux densities in the pointing 
where the sources are closest to the field center. The total flux density of each source in 
the mosaic and in the individual pointings are in very good agreement (the median value of 
the ratio is 1.03 with an r.m.s. dispersion of 0.04) confirming that the total flux density is 
properly recovered. On the other hand, as expected for observations affected by bandwidth 
smearing, the peak flux densities are underestimated in the final mosaic with respect to the 
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peak fluxes in the individual pointings where the sources are closest to the field center. To 
quantify this effect we selected only sources which are within a radius -Rmin < 5' from the 
field center in, at least, one individual pointing. The peak flux density of a source within 5' 
from the fleld center is almost unaffected (less than 5% for the VLA-COSMOS observations) 
by bandwidth smearing and can be compared with the peak flux in the mosaic. The image 
of the source in the mosaic is the combination of all the images in the individual pointings 
where the source is within the cut-off radius (16.8', see S07). 

In Figure [1] we show the peak flux ratios for compact sources (angular sizes less than 
4") with peak fluxes brighter than 0.2 mJy/beam and -Rmin < 5' as a function of the radial 
distance from the center in the flnal mosaic. 

The ratio is about 0.8 ± 0.05 for sources in the inner region of the mosaic, within 30' 
from the mosaic center. In the region between 30' and 45', the ratio is 0.9 ± 0.03, implying 
a smaller correction in this area. This is due to the fact that the number of overlapping 
pointings is smaller at larger radii ( see Fig. [2]). Within these two regions the effect appears 
to be independent from the position in the mosaic considering the dispersion of the points. 
Therefore, we adopted a constant bandwidth smearing correction factor of 0.8 for the sources 
within 30' from the mosaic center and a constant value of 0.9 for those between 30' and 45'. 

We performed this analysis in the inner 1 deg^ (1x1 deg. square region) of the flnal 
mosaic since this region has the deepest and most uniform sensitivity and it will be used to 
derive the source counts in the following Section. Beyond 45' from the mosaic center this 
correction is no more valid as each source in the mosaic receives a contribution from only one 
pointing and therefore the bandwidth smearing correction will be a function of the distance 
from the center of the individual pointing where each source is located, rather than of the 
distance from the mosaic center. 

Consequently, to derive the radio source counts in the inner 1 deg^ we divided all the 
peak flux densities by 0.8 or 0.9 depending on the distance from the center of the mosaic 
and derived the new resolved/unresolved criterion following the procedure described in S07. 
Since the ratio between total and peak fluxes is a direct measure of the extent of a radio 
source, we used it to discriminate between resolved or extended sources (i.e. larger than the 
beam) and unresolved sources. 

Figure [3] shows the ratio of the total integrated flux density ^totai and the corrected peak 
flux density S'peak as function of the signal to noise ratio. To select the resolved sources we 
have fltted a lower envelope in Fig. [3] which contains ~ 95% of the sources with S'totai < "S'peak 
and mirrored it above the S'totai /•Speak = 1 line. The criterion of 95% of sources contained in 
the lower envelope of Fig. [3] is slightly different from that used in S07, where a value of 99% 
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was adopted. We reckon that the value adopted here is the best compromise to properly 
classify the sources as resolved or unresolved because it fully characterises the shape of the 
distribution of sources with Stotai < •S'peak and, at the same time, it excludes the most extreme 
values. This envelope is described by the equation 

^total/^pcak = 1 + [l00/(^peak/rms)2-^] (1) 

It is important to note that in Fig.Othere is no systematic offset from the 5'totai/ 'S'pcak = 1 
line at high SNRs. This is an a posteriori confirmation that we adopted an appropriate 
bandwith smearing correction. Considering only sources with SNR> 5 in the inner 1 deg^, 
we obtain 484 resolved sources, for which the total flux is given by the total flux of the 
Gaussian fitting, and 1208 unresolved ones, for which the total flux is set equal to the peak 
flux. 

It is worth noting a few points regarding the effects of the correction on the VLA- 
COSMOS Large Project Catalog: 

• The correction has been determined only out to a radius of 45 arcmin from the mosaic 
center and should be safely adopted only in this region. 

• The catalog is still selected on the old, pre-bandwidth smearing correction, peak flux 
densities or signal-to-noise ratio, and we applied the bandwidth smearing correction 
only to the sources already in the catalog. In the remaining we will refer to the signal- 
to-noise ratio as the original one, calculated with the uncorrected peak flux density. 

• The correction affects the resolved/unresolved classification of a source and therefore 
also the total flux density and size for the sources, previously classified as resolved, 
that are found unresolved after correcting for the smearing. This affects about 500 
sources, almost one third of the total number of sources in the inner 1 deg^ with SNR 
> 5. 



3. Survey completeness 

3.1. The visibility area 

The VLA-COSMOS Large Project Catalog contains ~ 3600 radio sources extracted 
from a 2 deg^ area. As shown in S07, the noise is not uniform over the whole area and 
increases in the outer regions. The central 1 deg^ area has the lowest and most uniform 
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noise (mean rms 10.5 yuJy) and therefore it is very well-suited to derive the source counts. 
Extending the analysis to the whole area and catalog would increase the number of objects 
at higher flux densities (> 0.1 — 0.2 mJy), but it would not add much to the accuracy of 
the source counts at the lowest flux densities, which is the main focus of this paper. On the 
contrary, the correction factors that would need to be applied to take into account the loss of 
the faintest sources in the outer and noiser region would add more systematic uncertainties 
at the faint end of the source counts. Furthermore, in our following analysis we consider 
only sources with signal to noise ratio greater than or equal to 5 (the VLA-COSMOS Large 
Project Catalog has been selected with a signal to noise ratio threshold of 4.5) and without a 
flag for possible sidelobe residuals, in order to minimize contamination from spurious sources. 
Given these constraints, we have a catalog of 1692 sources with SNR > 5.0 within the inner 
1 deg^ region. 

The visibility area of the inner 1 deg^ of the VLA-COSMOS survey as a function of 
the rms noise is shown in Fig. HI This is the area over which a source with given peak flux 
density (5 times the noise in our case) can be detected. Since the noise in this region is very 
smooth the visibility area increases very rapidly and becomes fiat above an r.m.s. noise of 
~ 13 /iJy/beam. The same curve, computed over the entire 2 deg^, becomes flat only above 
an r.m.s. noise of ~ 35 /iJy/beam (see Fig. 13 in S07). This plot has been derived using the 
original peak flux densities and not those corrected for the bandwidth smearing since the 
catalog is still selected on the basis of the uncorrected peak flux. 

3.2. The resolution bias and the intrinsic angular size distribution of sub-mJy 

radio sources 

Besides the visibility area of each source, other correction factors need to be considered 
to estimate the completeness of the radio catalog and to derive the intrinsic source counts: 
the errors introduced by the fitting routines, the noise bias, and the resolution bias can all 
be modelled using simulated samples of radio sources. In particular, the resolution bias can 
be quite important since the beam of the VLA-COSMOS observations is 1.5" x 1.4" (e. g. 
about 15 times less in area than for similar surveys carried out at the same frequency). Since 
the completeness of the radio catalog is defined in terms of the peak fiux, while the source 
counts are derived as a function of the total fiux density, corrections need to be applied to 
take into account resolved sources with peak fluxes below the catalog threshold and total 
fiux densities above the nominal limit of the survey. 

Our strategy to estimate the correction factor has been to simulate samples of radio 
sources, with a realistic fiux density and angular size distributions, to insert these sources 



- 7- 



in the mosaic image and to recover these mock sources using the same technique used to 
produce the catalog. From the comparison between the number of sources added to the 
image and those effectively recovered we can estimate the correction factors for different flux 
density bins. The correct choice of the flux density distribution is rather straightforward, 
as a broken power law derived from the observed distribution is a reasonable assumption, 
consistent with the extrapolation from shallower surveys. The choice of the angular size 
distribution for the simulated radio sources is more uncertain. While there is good evidence 
that the median of the intrinsic source size at 1 mJy is about 2" and decreases with decreasing 
flux density at sub-mJy level, it is not known how exactly this happens and what is the best 
analytical relationship to model this behaviour. As we will show later, this is an important 
parameter in deriving the completeness of the VLA-COSMOS catalog. In the past different 
scaling relationships between the median ang ular size, ^msH, and the total flux d ensity, .S*, 
have been used, e.g. 6'med oc S"^ with m = 0.3 (IWindhorst et al.lll990l ) or m = 0.5 ([Richards 
2000h . 



More recently, new su rveys have provided so me results on the angular size distributions 
of the /iJy radio sources. iFomalont et al.l (l2006l ) obtained a complete sample of 289 radio 
sources with flux density > 40 /zJy and found that 64% of objects are unresolved with sizes 
less than 1.2 arcsec. MERLIN observations at higher angular resolution (~ 0.2") provided 
consistent results for a sample of 92 r adio sources, 80 of which with flux density between 40 
yuJy and 200 /iJy (IMuxlow et al.ll2005l ). The angular size distribution derived from the latter 
sample of objects, which are all resolved, shows a dominant narrow gaussian component 
centered at ~ 0".7 and a tail of sources (~ 20% of the sample) with angular sizes between 2" 
and 4". We derived the median of the angular size distribution for the sources observed by 
Muxlow et al. (2005) in the flux density range 40-90 /zJy. The 55 sources in this flux density 
range have a median angular size of 0".7 (with a r .m.s. dispersion of 0".4) and a median flux 
density of 54 /iJy. For comparison, iBondi et al.l (120031 ) found that the median angular size 
for sources in the flux density range 0.4-1.0 mJy is 1".8 with a median flux density of 0.56 
mJy. The value of m consistent with these median angular sizes at different flux densities 
is m = 0.4. It is worth noting that this value for m has to be considered as a lower limit 
since the angular size s given bv IMuxIow et al. (2005) are the largest angular sizes while the 
angular sizes given in Bondi et al.l ( 2003 ) are the FWHM. 



These results provide strong evidences that radio sources with flux density around 100 
/iJy have typically sub-arcsecond sizes even if the quahty of the data (e.g. large number of 
upper limits in the angular sizes derived in Bondi et al. 2003) and the statistics (the above 
result on m is derived from less than 100 sources) do not allow one to derive a flrm estimate 
of m. 
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Given these uncertainties, we decided to use a complementary method to derive the 
intrinsic angular size distribution, simulating mock samples of radio sources with different 
angular size distributions following a general power law 6 cx 5™ with the exponent m ranging 
from 0.2 to 0.5, and normali z ed to the integral angular size distribution derived from the 



VVDS survey in lBondi et all (l2003f ). 



We constructed a simulated sample of radio sources containing about 3100 radio sources 
down to a total flux density level of 30 /iJy. This allows us to count also sources with flux 
density below the limit which, because of positive noise fluctuations, might have a measured 
peak flux density above the threshold. Then, all the simulated sources were randomly in- 
jected in the inner 1 deg^ of the VLA-COSMOS mosaic and subsequently recovered using 
the same procedures adopted for the real sources and binned in flux density intervals. This 
was repeated for 3 different samples, and for four different angular size distributions (with 
m equal to 0.2, 0.3, 0.4, and 0.5) yielding a total of 12 simulated samples with more than 
35,000 simulated radio sources. Finally, from the comparison between the number of simu- 
lated sources detected in each bin and the number of sources in the simulated input sample 
in the same flux density bin we derived the correction factors for each value of the exponent 
m. These correction factors, listed in Table [T], account for the visibility area, noise bias, 
fitting errors and resolution bias. 

From Tabled] it is evident that for the brightest sources the correction factors are < 10% 
and the differences among the correction factors for the various angular size distributions are 
negligible. On the other hand, below 0.20 mJy the correction factors become quite large with 
significant differences for the various angular size distributions. Therefore, it is important 
to infer the intrinsic angular size distribution in order to reconstruct the properly corrected 
source counts. Another interesting feature of the correction factors which is worth noting is 
that the maximum is not at the lowest flux density bin but around 0.12 mJy. This is due to 
a combination of the resolution bias and the capability to classify a radio source as resolved, 
which depends on the signal-to-noise ratio. A source with a given total flux density and size 
will be catalogued in one of the following three classes: 1) undetected, if the peak flux is 
below the sample threshold, 2) detected with the correct total flux density if the signal to 
noise ratio is high enough to allow for proper deconvolution, or, 3) detected with a lower total 
flux density if the signal to noise ratio is not high enough to recognize the source as resolved. 
The sources in the last class are assigned a total flux density lower than the intrinsic one and 
equal to the peak flux. A significant fraction of the sources with intrinsic total flux densities 
< 0.15 mJy is redistributed to lower total flux densities since the signal-to-noise ratio is not 
sufficient to allow a proper deconvolution. 

To constrain the intrinsic angular size distribution we performed the following test. We 
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derived the median major axis of the simulated samples of radio sources for each of the tested 
power laws (m = 0.2 — 0.5) in the flux density range 0.25-0.4 mJy. This flux density range 
was chosen as a compromise to minimise the number of unresolved sources (and therefore 
upper limits on source size) and allow for a statistically significant large sample of sources 
reaching as close as possible to the lowest fluxes. We calculated both the median major axis 
of the original injected sources (row 1 in Tab. [2]) and the median major axis of the simulated 
sources recovered by the fitting routine (row 2 in Tab. [2]). These latter values were obtained 
considering also the upper limits derived for the sources fitted as unresolved (< 15% of the 
total number of sources in this flux density range). 

Table [2] shows that, in the flux density range 0.25-0.4 mJy, although the measured 
median major axis of the simulated sources is somewhat larger than that actually injected, 
the differences are not statistically significant for any given value of m. On the other hand, 
different values of m yield to significantly different median major axis. 

Therefore, the value of the median major axis in this flux density range for the real 
sources could be used to roughly discriminate between different intrinsic angular size distri- 
butions. Unfortunately, the measured sizes of the radio sources are affected by bandwidth 
smearing and they can not be used for such a test. However, as we have seen in Section 
2, the total flux densities are properly determined and the ratio between the total and the 
corrected peak flux densities can be used as an estimate of the area covered by the radio 
sources. Then we made the assumption that radio sources are extended in only one direction 
and we derived the major axis from the ratio between the total and the corrected peak flux 
densities. This assumption can be justified since extended radio emission is generally asso- 
ciated to elongated features (e.g. radio jets in AGN or spiral arms in star-forming galaxies), 
and when this is not true the derived source sizes will be an upper limit of the real ones. 
In this way we obtained a median value of the major axis for the observed sources in the 
flux density range 0.25-0.4 mJy of 0".99 ± 0".09. From the comparison with the values listed 
in Tab. [2l we conclude that the observed median major axis is more consistent with that 
obtained using simulated samples of radio sources following an angular size distribution with 
m = 0.5 or m = 0.4. 

With bot h meth ods, the comparison be tween the angular size distribution obtained by 



Muxlow et al.l (120051 ) and iBondi et al.l (120031 ) in different flux density intervals and the com- 
parison between the estimate of the true angular sizes and those derived from the simulated 
samples, we obtained consistent results which exclude values of m = 0.2 and m = 0.3 and 
suggest m ~ 0.5. In the following analysis we apply the correction factors obtained for 
m = 0.5 to the radio source counts. 
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The radio source counts 



We have derived the radio source counts at 1.4 GHz in the inner 1 deg^ region of 
the COSMOS field from the VLA-COSMOS catalog, corrected for the effect of bandwidth 
smearing and with a signal-to-noise threshold of 5. The source counts up to ~ 1 mJy are 
presented in Table [31 For each flux density bin we give the mean flux density, the observed 
number of sources, the corresponding differential source density dN/dS (in sr~^Jy~^), the 
normalised differential counts (dA^/dS')S'^'^ (in sr~^Jy^'^) with their estimated Poissonian 
errors, the incompleteness correction factors from Table [1] for m = 0.5, and the cumulative 
number of effective sources after applying the corrections. The correction factors above 0.2 
mJy are < 10%. Such variations are roughly consistent with the poissonian error of the 
source counts and for this reason the source counts above 0.2 mJy are assigned a correction 
factor of 1. The mean flux density in each bin was calculated as the geometric mean of the 
flux density extrema. 

In Fig. and Fig. [H] we show the raw and the corrected source counts derived from the 
VLA-COSMOS sur vey compared to thos e from other survey s at 1.4 GHz: the VL A-VVDS 
jBondi et all boosh . the ATCA-HDF_S diuynh et al.l \2003i]. t he VLA-HDF_N jRichards 
2000|), the PDF su rvev fEopkins et al.l l2003h. the SSA 13 held JPomalont et~aDl2006h and 
the FIRST survey ( IWhite et al.l 119971 ). The full range of source counts up to ^ 100 mJy is 
shown in Fig. O and a blow-up of the sub-mJy region is shown in FiglHl In both figures the 
raw source counts from the VLA-COSMOS survey are shown with empty circles and those 
corrected for incompleteness with filled ones. The solid line in Fig. [5] is a linear least-squares 
sixth-order polynomial fit obtained using the VLA-COSMOS source counts above 0.06 mJy 
supplemented with the FIRST source counts above 2.5 mJy and up to 1 Jy. The resulting 
polynomial fit is given by 



log[{dN/dS)/iS-'-')] = ^aaiog(5/mJy)]^ 



(2) 



i=0 



with ao = 0.805, ai = 0.493, = 0.564, = -0.129, = -0.195, ag = 0.110, 
and ag = —0.017. The residuals from the polynomial fit have an rms of a bout 0.06 in the 
logarithm of the normalized counts. As discussed by iHopkins et al.l (120031 ) the sixth-order 
polynomial is necessary to follow the different curvatures shown by the source counts. In 
Fi g. [5] we also show for c omparison with a dashed line the sixth-order polynomial fit obtained 
by iHopkins et al.l (120031 ) from the Phoenix Deep Survey source counts. The two fits are in 
reasonably good agreement with our fit slightly higher for fiux density below 0.2 mJy. 



It is evident that in some cases, the field-to-field differences in the radio source counts 
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below 0.5 mJy are larger than the combined errors. It is likely that these differences (at 
least partly) can be ascribed to instrumental defects or different recipes to account for the 
incompleteness. In particular, the assumption of the intrinsic angular size distribution of 
the fiJj radio sources in high resolution radio surveys is certainly a factor, as we have shown 
in the previous Section. Nevertheless, the most extreme cases (e.g. t he SSA 13 and HDF 
North fields) should reflect real cosmic variance (IFomalont et al.ll2006l ). 



There is a particularly good agreement between the source counts derived from the 
VLA-COSMOS, after correction, and the HDF-S fields. This is reassuring for the resolution 
bias correction we applied since the HDF-S field was observed with a resolution of about 5", 
and for this reason is much less affected by the precise shape of the intrinsic angular size 
distribution. 

It is worth noting that while the VLA-COSMOS survey is not as deep as the SSA13 and 
HDF-N surveys, it is, by far, the one with the most robust statistics. At a flux density level 
of ~ 65 /iJy the VLA-COSMOS survey counts almost 400 sources while typical numbers for 
other surveys, in this flux density range, are between 50 and 90 sources. 

A possible drop off in the radio source counts below ~ 100 /iJy has been already detected 
and discussed by other surveys imaging various fields with differe nt sensitivity and resolution 



(iHopkins et al.l l2003l : iHuynh et al.l l2005l : iFomalont et al.l l2006l . e.g.). Although, in every 



case, the drop off relies on the faintest bins it should not be immediately discounted as 
incompleteness in the source counts as it is accepted that th e flat region of the nor malised 
source counts below ~ 0.3 mJy can not continue indefinetely (jWindhorst et al.lll993f ). From 
the VLA-COSMOS radio source counts there is some indication that the density of sources is 
beginning to decrease below 100-150 fiJj producing a downturn in the observed radio source 
counts, even if the result is strongly dependent on the fidelity of the first flux density bin. 
It is important to note that if incompleteness is responsible for the drop off observed in the 
first flux density bin of Fig. [6l it would mean that we are missing about 130 sources in the 
range 0.06-0.0735 mJy, which is extremely unlikely. So, while incompleteness at the lowest 
flux density levels can still be present at some level, it can not explain by itself the drop off 
observed in the radio source counts. Thus, this decline can provide important limits when 
modeling the populations contributing to the radio source counts. 

As we have explained in the previous section, the intrinsic angular size distribution 
adopted to derive the incompleteness correction factor has a large impact in this flux density 
range, but adopting distributions with m < 0.5 would only enhance the drop off below ~ 0.15 
mJy. 



Also the uncertainties in the bandwidth smearing correction could, in principle, sig- 
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nificantly modify the radio source counts and the observed drop off observed at lower flux 
density. We have tested this assuming two extreme constant correction over the 

whole area of 0.8 and 0.9 respectively. As we can expect, the first case would produce a very 
small difference, well within the Poissonian error boxes, since we already applied the 0.8 
correction over almost 80% of the 1 deg^. On the other hand, adopting a correction factor 
of 0.9 over the whole area, would significantly lower only the radio source counts in the first 
three bins, increasing the slope of the drop off region. As we have shown in Section 2, the 
assumption of a constant correction over the whole 1 deg^ is not corroborated by the real 
data. The data require a mixture of the two correction factors and with the previous test we 
have just verified that the drop off in the source counts below 0.15 mJy can not be produced 
by the applied smearing corrections. 

We can conclude that given the large statistics and the fact that the incompleteness 
factors have been carefully investigated and applied, the decrease in the source counts below 
~ 0.10 — 0.15 mJy observed in the VLA-COSMOS radio source counts should be considered 
real. 



5. Summary 8z conclusions 

In this paper we have presented a further analysis of the VLA-COSMOS Large Project 
Catalog with the goals of: 1) quantifying and correcting for the bandwidth smearing effect 
which affects the VLA-COSMOS radio mosaic and catalog; 2) determining the incomplete- 
ness of the VLA-COSMOS radio catalog, through extensive use of Monte-Carlo simulations 
taking into account different intrinsic angular size distributions, in order to quantify the 
effects of noise bias and resolution bias; 3) deriving the radio source counts at 1.4 GHz in 
the VLA-COSMOS field. 

The bandwidth smearing correction factors have been determined from the comparison 
of peak and total fiux densities in the final mosaic and in individual pointings where the 
source is within 5' from the center. We have found that a two-value correction factor is a 
good approximation: out to a radius of 30' from the center, the peak flux densities in the 
VLA-COSMOS catalog must be divided by a factor 0.8, while in the region between 30' 
and 45' from the center by a factor 0.9. Inside these two regions the correction factor is 
rather constant, within the dispersion of the data, with a sharp transition from 0.8 to 0.9 
at a radial distance of 30' reflecting the change in overlapping pointings. Consequently, a 
new classiflcation in resolved and unresolved sources has been derived. A number of sources, 
originally classified as resolved, turned out to be unresolved after the bandwidth smearing 
correction was applied. For these sources also the total flux density was corrected. 
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We simulated several samples of radio sources with a flux density distribution compat- 
ible with the observed one and with different angular size distributions. The angular size 
distribution was modelled using 9 oc S"^ with m = 0.2, 0.3, 0.4, 0.5. These mock sources 
were added to the radio mosaic and recoverd using the same procedure that yielded to the 
radio catalog. In such a way wc were able to quantify the effects of the combined noise and 
resolution bias affecting the completeness of the radio catalog. A comparison of the median 
angular size at a flux density level not significantly affected by the resolution bias for the 
simulated samples of sources and the observed one, strongly suggest a general distribution 
of radio source sizes following 9 oc S^'^. 

Finally, we derived the radio source counts in the VLA-COSMOS field (both uncorrected 
and corrected for the incompleteness). The source counts extend down to 60 /xJy and are 
consistent with those derived by other surveys but with much more robust statistics. In 
particular, the drop off of the source counts below ~ 100 //Jy has to be considered real and 
not due to incompleteness. 

The National Radio Astronomy Observatory (NRAO) is operated by Associated Uni- 
versities, Inc., under cooperative agreement with the National Science Foundation. KJ 
acknowledges support by the German DFG under grant SCHI 536/3-1. 
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Table 1: Correction factors for different angular size distributions 



Flux bin 


m = 0.2 


m = 0.3 


m = 0.4 


m = 0.5 


(mJy) 










0.0600-0.0735 


1.29 


1.19 


1.06 


0.99 


0.0735-0.0900 


1.59 


1.35 


1.23 


1.16 


0.0900-0.1103 


2.12 


1.81 


1.53 


1.36 


0.1103-0.1351 


2.34 


2.12 


1.81 


1.60 


0.1351-0.1655 


1.60 


1.46 


1.32 


1.29 


0.1655-0.2028 


1.24 


1.21 


1.21 


1.15 


0.2028-0.2484 


1.05 


1.03 


1.00 


0.99 


0.2484-0.3043 


1.03 


1.02 


0.95 


0.93 


0.3043-0.4564 


1.10 


1.08 


1.09 


1.07 


0.4564-0.6846 


1.11 


1.10 


1.10 


1.09 


0.6846-1.0270 


0.98 


0.98 


0.97 


0.98 


> 1.0270 


1.07 


1.07 


1.08 


1.09 



Table 2: Median major axis in the flux density range 0.25-0.4 mJy for the input and output 
simulated samples 



Sample 


m = 0.2 


m = 0.3 


m = 0.4 


m = 0.5 




^maj,med 


^maj,med 


(^maj,med 


(^maj,med 


Sim. Input 


1.32 ±0.05 


1.17 ±0.04 


1.03 ±0.04 


0.91 ±0.04 


Sim. Output 


1.36 ±0.06 


1.24 ±0.06 


1.11 ±0.05 


1.00 ±0.05 
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Table 3: The 1.4 GHz radio source counts 



s 


<S> 


N 


dN/dS 


{dN/dS)S^-^ 


C 


Nc(> S) 


mJy 


mJy 




sr^-'^Jy"-'^ 


sr~^Jy^-^ 






deg-2 


0.0600-0.0735 


0.066 


380 


9.24 X 10^° 


3.32 ±0.17 


0.99 


1669 ±41 


0.0735-0.0900 


0.081 


324 


6.43 X 10^° 


3.84 ±0.21 


1 


.16 


1515 ±39 


0.0900-0.1103 


0.100 


203 


3.29 X 10^° 


3.26 ±0.23 


1 


.36 


1336 ±37 


0.1103-0.1351 


0.122 


142 


1.88 X 10^° 


3.09 ±0.26 


1 


.60 


1246 ± 35 


0.1351-0.1655 


0.150 


135 


1.46 X 10^° 


3.99 ±0.34 


1 


.29 


822 ± 29 


0.1655-0.2028 


0.183 


109 


9.61 X 10^ 


4.36 ±0.42 


1 


.15 


577 ± 24 


0.2028-0.2484 


0.224 


74 


5.33 X lO'^ 


4.02 ± 0.47 


1 


.00 


393 ± 20 


0.2484-0.3043 


0.275 


72 


4.23 X 10^ 


5.30 ±0.63 


1 


.00 


319 ±18 


0.3043-0.3728 


0.337 


45 


2.16 X lO'^ 


4.49 ± 0.67 


1 


.00 


247 ± 16 


0.3728-0.4566 


0.413 


33 


1.29 X lO'^ 


4.47 ± 0.78 


1 


.00 


202 ± 14 


0.4566-0.5593 


0.505 


28 


8.95 X 10^ 


5.14 ±0.97 


1 


.00 


169 ±13 


0.5593-0.6851 


0.619 


24 


6.26 X 10^ 


5.97 ± 1.28 


1 


.00 


141 ± 12 


0.6851-0.8393 


0.758 


12 


2.56 X 10^ 


4.05 ± 1.17 


1 


.00 


117± 11 


0.8393-1.0282 


0.929 


12 


2.09 X 10^ 


5.49 ± 1.58 


1 


.00 


105 ± 10 
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Fig. 1. — Ratio between the peak flux densities in the final mosaic and in the individual 
pointing where the source is within 5' from the center versus the radial distance from the 
center in the final mosaic. Only sources with peak flux density greater than 0.2 mJy/beam 
and sizes smaller than 4" are plotted. The straight lines at S (mosaic) /S (single pointing)= 0.8 
and 0.9 for R < 30' and R > 30' show the adopted bandwidth smearing correction factors. 
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Fig. 2. — Layout of the 23 pointings for the VLA-COSMOS observations. The two circles 
have a radius of 30' and 45'. 
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Fig. 3. — Ratio of the total flux to the corrected peak flux (see text for details) as a function 
of the signal-to- noise ratio. The solid hues show the upper and lower envelopes of the flux 
ratio distribution containing the sources considered unresolved. Open symbols show the 
sources considered extended. This Figure should be compared to Fig. 15 in S07. 
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Fig. 4. — Visibility area of the inner 1 deg^ of the COSMOS survey. 
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Fig. 5. — Radio source counts at 1.4 GHz from the VLA-COSMOS survey (dots) and from 
other surveys. Empty circles show the radio counts not corrected for incompleteness, filled 
circles the corrected ones using m = 0.5. The VLA-COSMOS source counts are shown along 
with those obtained by other deep surveys (see text). The solid line is least-squares sixth- 
order polynomial fit obtained usi ng the VLA-COSMO S and the FIRST source counts. The 
dashed line is the fit obtained by lHopkins et al.l (120031 ). 
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Fig. 6. — Blow-up of the radio source counts in the sub-mJy region. Symbols are the same 
as in Fig. O 



